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Quantum aspects, such as electron tunneling between closely separated metallic nanoparticles, 
are crucial for understanding the plasmonic response of nanoscale systems. We explore quantum 
effects on the response of the conductively coupled metallic nanoparticle dimer. This is realized 
by stretching a nanorod, which leads to the formation of a narrowing atomic contact between the 
two nanorod ends. Based on first-principles time-dependent density-functional-theory calculations, 
we find a discontinuous evolution of the plasmonic response as the nanorod is stretched. This is 
especially pronounced for the intensity of the main charge-transfer plasmon mode. We show the 
correlation between the observed discontinuities and the discrete nature of the conduction channels 
supported by the formed atomic-sized junction. 
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Metallic nanoparticles are of great interest be¬ 
cause they have potential in many applications, such 
as biosensor development [1], spectroscopy [2], and 
nanobiomedicine [3]. This is thanks to their versatile 
optical properties that are mainly determined by the 
localized surface plasmon resonances supported by the 
nanoparticle. These excitations are characterized by the 
collective oscillations of the delocalized electrons in re¬ 
sponse to light. Understanding the physics determining 
their dependence on nanoparticle parameters is thus a 
major subject of research seeking for different routes to 
tailor materials with desired optical properties. Along 
with the size, shape, environment, and atomic compo¬ 
sition of the individual nanoparticles, the plasmonic re¬ 
sponse can be modified by forming extended systems, 
such as nanoparticle arrays and agglomerates [4]. Among 
these, nanoparticle dimers are not only the simplest con¬ 
figurations for studying the plasmonic interparticle inter¬ 
action, they also serve as basic components in the design 
of nanoantennas [5, 6]. 

The plasmonic response of large nanoparticle dimers 
can be mostly understood by classical electrodynamics 
and the plasmon hybridization models [7, 8]. However, 
the classical description fails at subnanometer interpar¬ 
ticle separations, where quantum tunneling of electrons 
strongly affects the plasmonic response. Recently, this 
quantum regime has become experimentally accessible 
[9, 10] confirming several theoretical predictions [11-13]. 
Follow-up studies have further addressed quantum effects 
from dimers of mesoscopic sizes [14] down to the atom¬ 
istic details [15], also considering the formation of con¬ 
ductive contacts [16, 17], as well as molecular and atomic 
junctions between the nanoparticles [18-21]. 

In this Letter we report, based on first-principles cal¬ 
culations, that the plasmonic response of nanoparticle 
dimers shows important quantum features within the 
metal-contact regime. We consider a metallic nanorod 


subject to a stretching process leading to the formation 
of a narrowing atomic contact between the two nanorod 
ends until eventual breakage. Remarkably, the calcula¬ 
tions reveal a discontinuous evolution of the plasmon en¬ 
ergies and intensities during stretching. We find a strong 
correlation between the observed discontinuities and the 
quantized nature of the conductance of the narrow con¬ 
tact. This effect is pronounced for the main charge- 
transfer plasmon mode. 

For our study, we consider a flat-edged sodium nanorod 
obtained by extracting 25 atomic layers from the bee lat¬ 
tice along the (001) direction. The resulting system, 
shown in panel A of Fig. 1(a), comprises 261 atoms. The 
lattice constant was optimized to 4.14 A by energy min¬ 
imization of the nanorod (see below for the numerical 
details), leading to an approximate diameter of 13 A and 
an aspect ratio of 4. This size is close to the smallest 
metallic nanorods produced in recent experiments [22]. 
The use of sodium for modeling is not only computa¬ 
tionally convenient because of the single valence electron 
per atom needed in its description, it has also been suc¬ 
cessful in reproducing many general phenomena of sim¬ 
ple and noble metals in cluster science and plasmonics, 
both within atomistic [15-17, 23, 24] and continuum ap¬ 
proaches [11, 13, 25, 26]. 

The initial perfect nanorod is adiabatically stretched 
by pulling apart the frozen ends in steps of 0.05 A, as 
schematically indicated in panel A of Fig. 1(a). After 
each elongation step, the atomic positions of the sand¬ 
wiched central slab of 51 atoms are allowed to relax. The 
geometry relaxation is performed with density-functional 
theory [33] calculations using the GPAW code [34, 35]. We 
use the Perdew-Burke-Ernzerhof [36] (PBE) exchange- 
correlation (xc) functional. Selected snapshots of the 
simulated stretching process are shown in Eig. 1(a), la¬ 
beled from A to H. During the stretching, the break 
junction undergoes alternating sequences of elastic defor- 
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FIG. 1. Evolution of the plasmonic response of a sodium nanorod under stretching, (a) Selected snapshots of the simulated 
stretching process of the nanorod. Elongation distances d from A to H: 0, 6, 10, 14, 20, 22, 26, and 28 A. (b) Longitudinal 
photoabsorption spectra of the stretched nanorod as a function of energy hu. Solid lines following the plasmon modes (see 
main text for the label definitions) are drawn to guide the eye. (c) IID maps [see Eq. (2)] as a function of energy and z 
coordinate along the main axis of the nanorod. The maps have been smoothed along z to emphasize the gross features of the 
induced density [27]. Insets show the corresponding geometries (right) and spectra (top). Note that in panels E, E, and G, the 
low-energy part of the IID map and spectrum have been scaled by the indicated factor. 


mations and stochastic rearrangements. In this process, 
well-ordered structures can be formed in the junction, as¬ 
sociated to nanowires with stable magic radii [25, 37], a 
clear example of which is observed in panel F. Our sim¬ 
ulation is in line with previous approaches studying the 
breaking of metallic nanowires [24, 25, 38-40]. 

The optical properties of the described stretched 
nanorod are studied with time-dependent density- 
functional theory (TDDFT) [41] calculations. In brief, 
the time-dependent Kohn-Sham equations are solved in 
real time [42] with the adiabatic PBE xc functional as 
implemented in GPAW [43-45]. In this framework, the 
plasmons arise from the correlated time evolutions of 
the Kohn-Sham electrons in the mean field produced by 
other electrons. We consider only the longitudinal opti¬ 
cal response along the main axis (z) of the nanorod. The 
initial ground-state electronic structure of the system is 


perturbed by a weak delta-pulse electric field. The sub¬ 
sequent computation of the time evolution of the valence 
electrons yields the induced charge density Sp{r^t), from 
which the dynamic polarizability a{uj) can be derived. 
The response is characterized via the photoabsorption 
cross section of the system 

o-abs(w) = —Im[Q;(w)], (1) 

ceo 

where c is the speed of light and cq the vacuum permit¬ 
tivity. The plasmon resonances appearing in the pho¬ 
toabsorption spectrum can be further analyzed by the 
imaginary part of the integrated induced density (IID) 
defined as 


Sp{z^ uj) = Im 


/ 


Sp{r, uj) dxdy 


( 2 ) 


Note that the polarizability can be obtained from the IID 
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via Im[(a(co’)] oc f zSp(z,cj) dz. Thus, the IID serves as a 
2 ;-resolved spectrum summarizing the spatial distribution 
of the plasmon modes. 

We start our discussion by characterizing the different 
plasmon modes found in the calculated spectra shown in 
Fig. 1(b) for the increasing elongation distance d. The 
breaking of the narrowing junction after d = 26 A divides 
the spectra into two different regimes. We first focus 
on the contact regime (d < 26 A), characterized by the 
progressive quenching and emergence of several different 
plasmon modes. 

Initially (d = 0) the spectrum of the nanorod is gov¬ 
erned by a longitudinal plasmon (LP) mode at 1.5 eV. 
The dipolar character of this mode can be recognized in 
its corresponding IID in panel A of Fig. 1(c). Addition¬ 
ally, we observe a faint higher-order multipolar mode at 
2.9 eV (LP2). As the nanorod is stretched into a con¬ 
nected subnanorod dimer system, the evolved LP mode 
is better described as a charge-transfer plasmon (CTP), 
consistent with the literature on nanoparticle dimers [9- 
21, 46-48]. The former label, LP, emphasizes the single 
nanorod character of the mode. On the other hand, the 
latter label, CTP, reflects that the mode can only be sup¬ 
ported by the charge flow through the forming conduc¬ 
tive junction between the nanorod ends. Additionally, 
the LP2 mode of the single nanorod evolves into a bond¬ 
ing dipolar plasmon (BDP) in the dimer, both modes 
having a similar nodal structure in the IID. The BDP 
mode arises from the hybridization [8, 46] of the dipo¬ 
lar plasmon modes of the two subnanorods. In contrast 
to the CTP, the BDP mode exists via capacitive cou¬ 
pling, also when the junction is removed. For increas¬ 
ing d, the LP-CTP mode quenches accompanied by a 
broadening whereas the initially faint LP2-BDP mode 
grows, eventually becoming the dominant resonance in 
the response. Both peaks undergo strong redshifts, con¬ 
sistent with experiments and classical calculations per¬ 
formed on conductively connected metallic nanoparticles 
[48-51]. Thus, the redshifts can be ascribed to the com¬ 
bined effect of the simultaneous reduction of the diameter 
of the junction and the increasing length of the system. 

For d > 13 A, a second charge-transfer plasmon 
(CTP2) mode emerges at ^ 0.3 eV below the partially 
overlapping CTP peak. The CTP2 progressively grows 
with increasing d and eventually becomes comparable in 
strength to the CTP mode. A careful inspection of the 
IID maps in Fig. 1(c) (panels D-G) reveals subtle differ¬ 
ences between these two charge-transfer modes. While 
the CTP mode retains much of its dipolar character with 
a mostly localized IID in the far ends of the entire system, 
the CTP2 has a strong contribution from both the far and 
inner ends of the forming subnanorods. To further inves¬ 
tigate these two modes, we performed additional calcula¬ 
tions with different system sizes and junction geometries 
(not shown) suggesting that (i) the CTP2-like charge os¬ 
cillation is likely to be dominant for thin junctions and 


(ii) the two modes possibly merge in larger systems, in 
agreement with classical calculations where a single dom¬ 
inant charge-transfer mode is observed [49, 50]. 

Right before breaking, at d = 26 A, the stretched 
nanorod reaches the interesting limiting case of a 
monoatomic contact. In this configuration, the junc¬ 
tion is almost detached from the upper subnanorod [see 
panel G in Fig. 1(a)] and the CTP mode is strongly sup¬ 
pressed. In addition to the discussed charge-transfer plas¬ 
mon modes, a broad structure of comparable strength is 
observed at 1.6 eV. The corresponding IID map (G) in 
Fig. 1(c) reveals the contribution of the junction to the 
mode. Thus, this atomic configuration can be viewed as 
a heterogeneous trimer. 

After the breaking of the junction (d > 26 A), the 
spectrum changes abruptly. We observe that the CTP 
modes disappear completely. This is in contrast with the 
progressive onset of the CTP modes supported by quan¬ 
tum tunneling when bringing two metallic nanoparticles 
into contact [9, 10]. In our case, the quantum tunneling 
is rapidly suppressed after the breaking because of (i) 
the contraction of the broken junction towards the sub¬ 
nanorods, thus increasing the vacuum gap, and (ii) the 
small available tunneling area in the resulting geometry. 
Another major change in the spectrum is the appearance 
of a second bonding dipolar mode (BDP2) at 1.3 eV. 
From the IID map shown in panel H of Fig. 1(c) we ob¬ 
serve that this mode is associated with a strong dipolar 
charge oscillation within the attached remnant junction 
and the lower subnanorod, inherited from the weak pre¬ 
breakage feature at 1.6 eV discussed above. 

Having characterized the plasmon modes in the sys¬ 
tem, we proceed to analyze in detail their evolution. For 
this aim, we show in Figs. 2(a) and 2(b) the extracted en¬ 
ergies and intensities of the plasmon modes as a function 
of the elongation [52]. These figures reveal the rich struc¬ 
ture in the evolution of the modes. The stepped evolution 
of the different plasmon modes is particularly remarkable. 
This is especially pronounced for the CTP mode that un¬ 
dergoes strong intensity drops during stretching at many 
elongation values. We also observe discontinuous shifts 
in the plasmon energies at the same elongations, most 
notably at d = 7.5 A where the BDP and CTP modes 
redshift by 0.2 eV and 0.1 eV, respectively. The disconti¬ 
nuities can be traced back to the rearrangements in the 
junction, but such strong intensity changes can hardly 
be explained by the geometrical change alone. Indeed, 
the discontinuous evolution of the intensity of the CTP 
mode is very similar to the quantum features shown by 
the conductance curves measured during the breaking of 
metallic nanowires [53]. 

To examine the connection between the observed dis¬ 
continuities and the electron transport properties of the 
junction, we have calculated the zero-bias conductance 
of the junction by the non-equilibrium Green’s function 
technique [54] as implemented in GPAW [55, 56]. The con- 
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FIG. 2. Detailed analysis of the plasmon modes during 
the nanorod stretching, (a) Peak energy and (b) integrated 
intensity (the area under the peak) of the plasmon modes as a 
function of the nanorod elongation d [52] . The intensities are 
normalized so that the full spectrum integrates to the number 
of valence electrons (261). The line labels correspond to the 
labeling used in Fig. 1(b) and labels A-H to the geometries 
shown in Fig. 1(a). (c) Zero-bias conductance of the junction 
formed during the stretching. Horizontal lines mark values of 
1, 3, and 6 Go. Vertical lines mark the discontinuities found 
in the evolution of the plasmon modes. 

ductance results plotted in Fig. 2(c) show the well-known 
features arising in breaking atomic-sized junctions from 
the quantization of the conduction channels [53]: Upon 
stretching, in addition to the gradual evolution of conduc¬ 
tance (reflecting changes in conduction channels’ trans¬ 
mission probabilities), clear jumps of integer multiples of 
the quantum conductance unit Go = 2e^//i take place 


(due to the closing of channels). When these steps are 
compared with the discontinuities found in the evolution 
of the plasmon modes, we observe a nearly one-to-one 
match. Thus, we conclude that upon the stretching of the 
nanorod, the number of conduction channels supported 
by the junction is reduced in a discrete manner [27] af¬ 
ter the atomic rearrangements in the junction; this leads 
to abrupt changes in the allowed charge flow between 
the subnanorods as shown by the conductance, which is 
further reflected in the plasmonic response. This phase- 
coherent transport effect is particularly pronounced in 
the CTP absorption intensity, as expected from the char¬ 
acter of the mode. 

While the intensity of the CTP mode closely follows 
the trends shown by the conductance of the junction, 
even reflecting the increase of the conductance at d ~ 
21 A (ascribed to enhanced transmission due to the for¬ 
mation of the well-ordered junction geometry), the CTP2 
mode behaves differently. This is attributed to the previ¬ 
ously described differences in the spatial characteristics 
of the two modes: The reduction of the effective cross sec¬ 
tion of the junction along with the elongation uncovers 
the inner surfaces of the subnanorods necessary to sup¬ 
port the CTP2 mode and leads to its growth. The BDP 
mode, also having a strong localization in the inner ends 
of the subnanorods, is similarly affected. Further fine de¬ 
tails in the evolution of the spectra can be attributed to 
relatively complex stochastic atomic configurations. 

In conclusion, we performed a first-principles study of 
the plasmonic response of a stretched metallic nanorod. 
We found a quantized evolution of the plasmon modes 
and established a connection between the plasmonic dis¬ 
continuities and the discrete number of conduction chan¬ 
nels supported by the break junction. Thus, this phe¬ 
nomenon can be considered the plasmonic counterpart of 
the quantization of the conductance in atomic-sized junc¬ 
tions. The effect is most prominent in the main charge- 
transfer plasmon of the considered sodium nanorod. In 
noble metal systems the damping due to the d-electrons 
is weak for such a low-energy plasmon mode [58, 59], and 
thus, we expect similar quantized evolution to be observ¬ 
able also in these materials. As our findings are based on 
the formation of a thin nanocontact, they are most likely 
to be experimentally reproducible in a contact breaking 
situation [51, 60] of nanoparticle dimers rather than in a 
jump to contact one [9, 10, 16, 17, 48]. 
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SI Smoothing of the integrated induced density maps 

For completeness, we show in Fig. SI the same data as in Fig. 1(c) of the Letter, but 
without applied smoothing along the ^ coordinate. The smoothed figure emphasizes 
the envelope of the induced density that is important for characterizing the plasmon 
modes. In contrast to the smoothed data, the raw data shows rapid oscillations 
along the 2 ; coordinate. These oscillations can be attributed to size effects and 
charge screening at the surfaces (Friedel-like oscillations) [1-4]. 
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Figure SI: IID maps with no smoothing along the z coordinate. 
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S2 Junction electronic structure and transport properties 

In this section we provide a more detailed analysis of the junction electronic struc¬ 
ture leading to the quantization effects observed in our study. To this end we 
calculate the junction projected density of states (DOS) and show the quantized 
electronic structure giving an estimation of the conduction channels contributing 
to the charge transfer plasmons found in the stretched nanorod. 

We define the junction projected DOS (jDOS) as the projection of the density 
of states onto the junction, 


9iunct{E) = '^WiS{E - €i), (SI) 


where Wi and are the weight and eigenenergy, respectively, of the ith Kohn-Sham 
electron state with the wave function 'ipi{r). The weight Wi is defined as [5] 


/ 

J Z^[zq — ^Z,Zq-\-^Z^ 


|■^/)j(r)pdr 



(S2) 


The first integral on the RHS of Eq. (S2) is calculated within a slab of width 2Az 
centered at zq, the position of the plane with the smallest junction cross section. 
Here, Az is chosen as 4 A and zq is obtained from the minimum of the x, y- 
integrated ground-state electron density n^^y{z) = f f n(r)dxdy, the atomic-scale 
variations of which have been smoothed along the 2 ; direction. 

In Fig. S2 we show the calculated jDOS for the geometry F of the stretched 
nanorod [see Fig. 1(a) in the Letter]. For visualizing the results we have substituted 
the Dirac delta function S{E — ei) in Eq. (SI) by a normalized Gaussian function 

G,{E - ei) = - 2 — ( 33 ) 

a\/27r 

In Fig. S2, the results are shown for smoothing parameter values a of 0.15 eV and 
0.03 eV. By considering first the jDOS with a = 0.15 eV, we observe two asymmetric 
peaks with long tails towards higher energies. Their shapes are characteristic for a 
one-dimensional electron gas (IDEG) with the DOS g{E) r\j E [7]. This hints 
that we could model the junction as a idealized cylindrical confining potential of a 
length L and a diameter D, where L is large compared to D. The single-electron 
states in such a potential are characterized by the quantum numbers (m, n) where 


Junction projected DOS of geometry F 
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Figure S2: The junction projected density of states of the geometry F (elongation 
22 A). The Fermi level is at 0 eV. Note the characteristic DOS shape of a one¬ 
dimensional electron gas. The thin bars indicate the weights Wi of the individual 
states. 
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m = 0, ±1, ±2,... is the angular quantum number and n = 1, 2, 3,... the radial 
quantum number. Indeed, we can confirm the validity of this idealized picture 
by visually inspecting the wave functions contributing most to the two labeled 
peaks in jDOS. In general, the wave functions extend over the whole stretched 
nanorod, but by focusing at the junction region, we can associate the lower peak 
( 1 ) to a state corresponding to (m,n) = ( 0 , 1 ) and the higher peak ( 2 ) to two 
nearly-degenerate states corresponding to (m,n) = (±1,1). With the increased 
energy resolution (<j = 0.03 eV) we can distinguish the discrete states yielding the 
characteristic DOS shape. By visually inspecting the wave functions of the states 
at the long tails of the characteristic IDEG DOS, we can observe that these states 
correspond to the same (m, n) numbers but have different numbers of nodes along 
2 ; and/or they are mixed with different states at nanorod ends. Thus, the present 
metallic junction coupled with the nanorod ends through metal bonding appears 
to be long enough to support (quantized) one-dimensional electron “bands” over its 
length. These “bands” serve as conduction channels for the phase-coherent charge 
transport through the junction. This situation differs from molecular junctions [ 8 , 
9] which might have very different electronic structures and coupling mechanisms 
to the nanoparticles forming the dimer. 

In Fig. S3 we show the evolution of the jDOS for the occupied states during the 
nanorod elongation. We observe clear maxima in the jDOS which are associated 
to states having stronger density at the junction. These states generally rise in 
energy as the junction shrinks during the elongation. By visually inspecting wave 
functions that contribute most to the jDOS peaks, we can associate the states to 
the following (m, n) states of the cylindrical symmetry: 

1 ) (m^n) = ( 0 , 1 ) (1 state), 

2 ) {m,n) = (± 1 , 1 ) (2 states), 

3) (m, n) = (±2,1) and (m, n) = (0, 2) (close in energy, 3 states in total), 

4) (m^n) = (±3,1) (2 states), and 

5) {m,n) = (±1,2) (2 states). 

Note that the states with m 7 ^ 0 are not exactly degenerate in the real non- 
uniform system. The energies of the states closely coincide with those of a uniform 
cylindrical jellium (Na) nanowire of similar diameter [10]. All in all, this statistical 
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Figure S3: The junction projected density of states as a function of the nanorod 
elongation [Eqs. (S1)-(S3) with a = 0.15 eVj. The Fermi level is at 0 eV. Note that 
the the DOS has been multiplied by a factor of two for elongations larger than 
12 A. The dotted black lines track the general evolution of the main peaks labeled 
1-5. The vertical gray lines and the geometry labels A-H are consistent with Fig. 2 
in the Letter. 
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picture shows good agreement with the idealized model presented above in spite of 
the variations in the exact junction geometry during the stretching. 

The JDOS together with the presented model provides direct estimate on the 
number of open conduction channels supported the junction. By comparing the 
numbers of “bands” to the calculated conductance [see Fig. 2(c) in the Letter or 
Fig. S4 below], we observe that the numbers serve as upper values for the conduc¬ 
tance in units of Gq. In the elongation region denoted by C in Fig. S3, the total 
number of channels is estimated to be 1 + 2 + 3 = 6; in the regions D-F 1 + 2 = 3, 
and finally, for the single atom contact in the region G, there is one single channel. 
The atomistic nature and non-uniformity of the junction affects the transmission 
probability of each conduction channel, and thus the conductance values appear 
generally below these upper limit estimates. 

Next, we compare the transport properties of the junction, determined by the 
number of conduction channels and their transmission probabilities, to the charge 
transferred through the junction at the CTP energy. To this end, we have estimated 
the induced charge at the CTP mode as 

Sq{iUcTp) = / ^p{r,uJcTp) dr, (S4) 

J 2;e[— oo,2;o] 

where dp{r,uJcTp) is the induced charge density at the CTP energy cjctp, and the 
integration limit zq is the position of the plane with the smallest junction cross 
section as discussed above. In other words, the integration is carried over one of 
the nanorod ends. In Fig. S4 we show the induced charge on top of the calculated 
zero-bias conductance [see Fig. 2(c) in the Letter]. The agreement between these 
two quantities is remarkable even though the frequency-independent zero-bias con¬ 
ductance is compared to the light-induced charge. The agreement is clearest after 
the first discontinuity, after which the thin junction becomes the restrictive factor 
determining the properties of the CTP mode. 

Finally, we briefly note that the induced charge is directly related to the optical 
linear response of the system discussed in the Letter. The imaginary part of the 
dynamic polarizability, calculated as the Fourier transform of the time-dependent 
dipole moment normalized to the perturbing electric field strength, can be esti¬ 
mated as Im[(a(ci;cTp)] oc Llm[Sq{uJcTp)]^ where L is the effective length of the 
dipole over which the charge is oscillating at the CTP energy. 



Elongation (A) 


Figure S4: Comparison between the conductance and the induced charge at the 
CTP mode energy. Black line with circles: zero-bias conductance of the junc¬ 
tion during the stretching [from Fig. 2(c) in the Letter]. Red line with squares: 
the absolute value of the induced charge \Sq{uJcTp)\ (~ Iin[<^^(<^CTp)]) during the 
stretching. The vertical and horizontal gray lines are consistent with Fig. 2(c) in 
the Letter. 
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